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Abstract 

Theories of anti-commuting scalar fields are non-unitary, but they are of interest both 
in statistical mechanics and in studies of the higher spin de Sitter/Conformal Field Theory 
correspondence. We consider an Sp{N) invariant theory of N anti-commuting scalars and 
one commuting scalar, which has cubic interactions and is renormalizable in 6 dimensions. 
For any even N we find an IR stable fixed point in 6 — e dimensions at imaginary values of 
coupling constants. Using calculations up to three loop order, we develop e expansions for 
several operator dimensions and for the sphere free energy F. The conjectured F-theorem 
is obeyed in spite of the non-unitarity of the theory. The 1/N expansion in the Sp{N) 
theory is related to that in the corresponding 0{N) symmetric theory by the change of sign 
of N. Our results point to the existence of interacting non-unitary 5-dimensional CFTs 
with Sp{N) symmetry, where operator dimensions are real. We conjecture that these CFTs 
are dual to the minimal higher spin theory in 6-dimensional de Sitter space with Neumann 
future boundary conditions on the scalar field. For N = 2 we show that the IR fixed point 
possesses an enhanced global symmetry given by the supergroup OSp{l\2). This suggests 
the existence of OSp{l\2) symmetric CFTs in dimensions smaller than 6. We show that the 
6 — e expansions of the scaling dimensions and sphere free energy in our OSp{l\2) model are 
the same as in the g —)■ 0 limit of the g-state Potts model. 
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1 Introduction and Summary 

Among the classic models of quantum field theory, a prominent role is played by the 0{N) 
invariant theories of N massless scalar fields 0*, which interact via the potential For 

any positive N these models possess interacting IR fixed points in dimensions 2 < d < 4 [^ . 
These theories contain 0{N) singlet current operators with all even spin, and when N is 
large the current anomalous dimensions are ~ 1/N [^. Since all the higher spin currents 
are nearly conserved, the 0{N) models possess a weakly broken higher spin symmetry. In 
the Anti-de Sitter/Conformal Field Theory (AdS/CFT) correspondence [3]-[^, each spin s 
conserved current in a d dimensional CFT is mapped to a massless spin s gauged field in 
d + 1 dimensional AdS space. For these reasons, it was conjectured that the singlet sector 
of the critical 0{N) models in d = 3 is dual to the interacting higher spin theory in d = 4 
containing massless gauge fields of all even positive spin [7 10 . This minimal Vasiliev theory 
also contains a scalar field with = —2./i\j^g, and the two admissible boundary conditions 


on this field 11distinguish the interacting 0{N) model from the free one (in the latter 
all the higher spin currents are conserved exactly). A review of the higher spin AdS/CFT 


dualities may be be found in 13 


A remarkable feature of the Vasiliev theories (7-10 is that they are consistent not only 
in Anti-de Sitter, but also in de Sitter space. On general grounds one expects a CFT dual to 
quantum gravity in dS '4 to be a non-unitary theory defined in three dimensional Euclidean 


space 14 . In 15 it was proposed that the CFT dual to the minimal higher spin theory in 


dS'4 is the theory of an even number N of anti-commuting scalar fields y* with the action 




( 1 . 1 ) 


This theory possesses Sp{N) symmetry, and Qij is the invariant symplectic matrix. This 
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model was originally introduced and studied in 16,17 , where it was shown to possess an 


IR hxed point in 4 — e dimensions. The beta function of this model is related to that of the 


0{N) model via the replacement N —>■ —N. According to the proposal of 15 , the free UV 


hxed point of (1.1) is dual to the minimal higher spin theory in dS 4 with the Neumann future 
boundary conditions on the = ‘2/i‘^g bulk scalar held, and its interacting IR hxed point to 
the same higher spin theory but with the Dirichlet boundary conditions on the scalar held. 
In the latter case, the higher spin symmetry is slightly broken at large N, and the de Sitter 
higher spin gauge helds are expected to acquire small masses through quantum ehects. A 
discussion of the de Sitter boundary conditions from the point of view of the wave function 
of the Universe was given in 


21 


22 


In this paper we consider an extension of the proposed higher spin dS/CFT correspon¬ 


dence 15 to higher dimensional de Sitter spaces, and in particular to dS^. Our construction 


mirrors our recent work 23,24 on the higher dimensional extensions of the higher spin 


AdS/CFT. It was observed long ago 25 -27 that in d > 4 the quartic 0{N) models possess 
UV hxed points which can be studied in the large N expansion. The UV completion of the 


0{N) scalar theory in 4 < d < 6 was proposed in 23 ; it is the cubic 0{N) symmetric theory 
of V -|- 1 scalar helds a and 0h 


5 = I ‘t‘x{ 


( 1 . 2 ) 


For sufficiently large N, this theory has an IR stable hxed point in 6 — e dimensions with 
real values of gi and g 2 - The beta functions and anomalous dimensions were calculated 
to three loop order [^,24 ,^and the available results agree nicely with the 1/N expansion 


of the quartic 0{N) model at its UV hxed point 29-33 when the quartic 0{N) model is 


continued to 6 — e dimensions. The conformal bootstrap approach to the higher dimensional 


0{N) model was explored in 34-36 


To extend the idea of 15 to dS^+i with d > 4, we may consider non-unitary CFTs (1.1) 


which in d > 4 possess UV hxed points for large N. The 1/N expansion of operator scaling 
dimensions may be developed using the generalized Hubbard-Stratonovich transformation, 
and one hnds that it is related to the 1/N expansion in the 0{N) models via the replacement 


potential difficulty with this picture is that unitarity in space requires that a massive field of 

spin s > 1 should satisfy m? > {s — l){s + d — 18 -20 . In other words, there is a finite gap between 


massive fields and massless ones (the latter are dual to exactly conserved currents in the CFT). However, 
since the masses are generated by quantum effects and are parametrically small at large N, this is perhaps 
not fatal for bulk unitarity. It would be interesting to clarify this further. 


■^The one loop beta functions of the model (1.2) in d = 6 were first calculated in 28 . 
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N —>■ —N. In (i = 5 this interacting fixed point should be dual to the higher spin theory in 
dSe with Neumann boundary conditions on the = Q/i'^s scalar field (corresponding to the 
conformal dimension A = 2 + 0{1/N) on the CFT side)|^ In search of the UV completion 
of these quartic CFTs in 4 < d < 6, we introduce the cubic theory of one commuting real 
scalar field a and N anti-commuting scalar fields yh 

S = J d^x(^nijd^x'd^X^ + ^ ^ • (1-3) 


Alternatively, we may combine the fields y® into N/2 complex anti-commuting scalars d", 


a = 1,..., N/2 [16, 17 ; then the action assumes the form 


S = j d'^x(^d^9<^d>^9° + ^ (d^af + gia9^9^ + ^g2(T^^ ■ (1-4) 

We study the beta functions for this theory in d = 6 — e and show that they are related to 


the beta functions of the theory (1.2) via the replacement N —>■ —N. For all N there exists 


an IR fixed point of the theory (1.3) with imaginary values of gi and g 2 - This is similar to 


the IR fixed point of the single scalar cubic field theory (corresponding to the N = 0 case of 


our models), which was used by Fisher 37 as an approach to the Lee-Yang edge singularity. 
The fact that the couplings are purely imaginary makes the integrand of the path integral 
oscillate rapidly at large a; this should be contrasted with real couplings giving a potential 
unbounded from below. 


Our results allow us to study the 6 — e expansion of the theory (1.3) with arbtrary N, 


and we observe that at finite N there are qualitative differences between the Sp{N) and 
0{N) models which are not seen in the 1/N expansion. In fact, for the Sp{N) model there 
is no analogue of the lower bound Ncrit that was found in the 0{N) case [^. For the 
lowest value, N = 2, we observe some special phenomena. In this theory, which contains 


two real anti-commuting scalars, it is impossible to formulate the quartic interaction (1.1); 


thus, the cubic lagrangian (3.1) seems to be the only possible description of the interacting 


theory with global Sp{2) symmetry. Furthermore, it becomes enhanced to the supergroup 
OSp{l\2) because at the IR fixed point the two coupling constants are related via g^ = 2gl. 
The enhanced symmetry implies that the scaling dimensions of a and y® are equal, and we 
check this to order e^. An example of theory with OSp{l\2) symmetry is provided by the 


^The d = 5 free Sp{N) model corresponds to Dirichlet scalar boundary conditions in the dSg dual. It 
should be possible to extend this free Sp{N) model/d^ higher spin duality to general dimensions, since the 


Vasiliev equations in {A)dSd+i are known for all d 10 . 
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g —0 limit of the g-state Potts model 38 We show that the 6 —e expansions of the scaling 
dimensions in onr OSp{l\2) symmetric theory are the same as in the g —)■ 0 limit of the g- 
state Potts modelj^ This provides strong evidence that the OSp{l\2) symmetric IR hxed 


point of the cnbic theory (3.1) describes the second order transitions in the ferromagnetic 


g = 0 Potts model, which exist in 2 < d < 6 39 . 

Using the resnlts of [40|, we also compnte pertnrbatively the sphere free energies of the 


models (1.3). The sphere free energy for the OSp{l\2) symmetric model is fonnd to be the 


same as for the g = 0 Potts model. In terms of the qnantity F = sin(^) log which 


was introdnced in 40 as a natnral way to generalize the F-theorem 411443 to continnous 


dimensions, we hnd that the RG flow in the cnbic Sp{N) models in d = 6 — e satishes 


Fuv > -^iR for all N > 2. We show that the same resnlt holds in the model (1.1) in d = 4 — e. 
0 This is somewhat snrprising, since for non-nnitary CFTs the ineqnality Fuv > -^iR is 
not always satished. It wonld be interesting to nnderstand if this is related to the “psendo- 


nnitary” strnctnre discnssed in 17 , and to the fact that these models are presnmably dnal 


to nnitary higher spin gravity theories in de Sitter space. 


2 The IR fixed points of the cubic Sp(N) theory 


The beta fnnctions and anomalons dimensions for the Sp{N) symmetric model (1.3) can 


be obtained by replacing N —)■ —N in the corresponding resnlts for the cnbic 0{N) model 


(1.2), which were compnted in 23,24 to three loop order. Indeed, writing the action in the 


complex basis (1.4), we see that the Feynman rnles and propagators are identical to those 


of the 0{N) theory written in the U{N/2) basis, the only difference being that the N/2 
complex scalars are anticommnting. Hence, for each closed loop of the d" we get an extra 
minns sign, thns explaining the replacement N —)■ —N. 


Using the resnlts in 23,24 , the beta fnnctions for the Sp{N) model are then fonnd to 


^We are grateful to Giorgio Parisi for informing us about this and for important discussions. 

®We thank Sergio Caracciolo for suggesting this comparison to us and for informing us about the paper 


39 


°The models (1.1) and (1.3) also satisfy the P"-theorem to leading order in the large N expansion (for all 
d), since the leading order correction to Ujv ~ Gr is of order and was shown to satisfy the F-theorem 
in the corresponding unitary 0{N) models 40 43 44 . 
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be: 


A = - 2 ^'! “ i 2 ( 47 r) 3 ^i + 125(15(2 - 52) 


1 


432(4 


:7r, 


.51 ((536 - 86iV)5t + 12(30 + lliV)53^2 + (628 - ^ 24515 ^ _ 13 ^ 4 ^ + 


P 2 = - ^92 + {^Ngl - Nglg^ - ?,gl) 


144(4 


:7r 


(2iN<j\ + 322Ngig^ + eONgfgl - SlNglgl - 125gl) +.... 


( 2 . 1 ) 

We have omitted the explicit three loop terms, which can be obtained from 23,24 . Similarly, 
the anomalous dimensions of the helds y* and a take the form 


“6(1^ + 432^L)6 “ ^^^ 2 ) + • • •, 


1 ?? 


la = 


Ngl - gj 1 

12(47r)3 432(471)6 


{2Ngl + 48iV5^52 - glgl - l?>gl) + 


( 2 . 2 ) 


With the beta functions at hand, we can look for non-trivial hxed points of the RG flow 
satisfying /9i(5i,52) = O 5 /S 2 (fi'i, 1 / 2 ) = 0- all positive iV, we hnd two physically equivalent 
hxed points with purely imaginary coupling constants, hence all operator dimensions remain 
real. These hxed points are IR stable for all N (the stability matrix Mij = ^ has positive 


eigenvalues). Note that this is diherent from the 0{N) versions of these models 23,24 


where one hnds a critical N, whose one loop value is ~ 1038, below which the IR stable hxed 
points with real coupling constants disappear. However, to all orders in the 1/N expansion, 
the hxed point couplings and conformal dimensions in the Sp{N) models are related to the 
ones in the 0{N) models by the replacement N —)■ —N. 

Figure shows the RG how directions for N = 2. The arrows indicate how the coupling 
constants how towards the IR. The two IR hxed points are physically equivalent because 
they are related by gi —)■ — 5 ^. At higher values of N, the qualitative behavior of the RG 
hows and hxed points remain the same. We still have a UV Gaussian hxed point, and two 
stable IR hxed points. 

A special structure emerges for N = 2. In this case we hnd the hxed point solution 


* _ o * * _ • / A ,67 2^ 

92 - 25i , 9i -^\| —^— ( 1 + + 0{e ^ 


(2.3) 


and the conformal dimensions of the fundamental helds are equal (the three loop term given 
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Figure 1: The zeroes of the one loop (3 functions and the RG flow directions for the 05^(112) 

model. The coordinates are dehned via gi = g 2 = and the red dots 

correspond to the stable IR hxed points. 


below can be obtained from the results in 23,24 


= 2 - —e-- 

^ 15 450 


269 - 702C(3) 
33750 


+ 


(2.4) 


We show in the next section that the equality of dimensions is a consequence of a symmetry 
enhancement from Sp{2) to the supergroup OSp{l\2). 

It is natural to ask whether symmetry enhancement can occur at other values of N. For 
instance, we can explicitly check for which N the dimensions of a and y are equal. A direct 
calculation using the beta functions and anomalous dimensions up to three loops shows that 
this only happens for N = 2 and N = —1. The latter case corresponds to the 3-state Potts 
model hxed point of the the theory with two commuting scalars (^: it has gl = —g 2 and 
enhanced symmetry]^ 

The anomalous dimensions of some composite operators may be similarly obtained from 


the results in 23,24 . Let us quote the explicit result for the quadratic operators arising from 


^For = — 1, there is also a (non-unitary) solution with gl = which corresponds to two decoupled 


Fisher models 37 , and hence the dimensions of the fundamental fields are trivially equal. 
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the mixture of and Vti-x^'X^■ These operators have the same classical dimension, so we 


expect them to mix. The 2x2 anomalous dimension mixing matrix 7 “^ was given in 23 ^4 
up to one-loop order. Extending those results to two loops, we hnd the mixing matrix 


gf(N+4) I 3?(fff(21Af-134)-6si52(2W+5)+55|) 

I 


6(47r)^ 108(47r)® 

gl(g2-6gi) _ gl(6g3(llAr+20)-gfg2(lljV-324)+12gig2_i3g3^ 
6(47r)^ 216(47r)® 


Ngi{6gi-g2) , (l8gi+161gig2+15g|) 

6(47r)^ ' 108(47r)® 

Ngf+4g^ A^gf(160gf+36gig2-41g|)-181gf 

6(47r)3 216(47r)6 


(2.5) 


where index ‘1’ corresponds to the operator Qijx'‘x\ ^ind index ‘2’ corresponds to The 
diagrams contributing to the calculation are listed in Fig. 

'OK *"00 * 'OK * 'KKCK* KOK" 'OK "OK" ■KKK(" 'KK( " 'OK 


Figure 2: Diagrams contributing to the mixing of cr^ and QijX^'X^ operators to two loop order. 


The two eigenvectors of 7 “^ give the two linear combinations of cx^ and VLijX^X^ 1 the 
eigenvalues x± give their anomalous dimensions, so that A± = d — 2 -|- 7 ±. We hnd that 
one of these combinations is a conformal primary, and the other one is a descendant of a. 
Indeed, after plugging in the hxed point couplings, we hnd A_ = A^. -|- 2. For instance, in 
the large N expansion we hnd the results 


A_ — Act -|- 2 — A -\- 


-40e + 6800e - 


N 


+ 


A_|_ — A -\- 


lOOe - 


N 


+ 


-49760e + 


237476 ^2 


-I- . . . , 


( 2 . 6 ) 


A2 


+ .... 


Upon sending N —)■ —N, the dimension A+ can be checked to be in agreement with the 
available large N results for the critical exponent u in the quartic 0{N) theorjj^ 32,45 . 


2.1 Estimating operator dimensions with Fade approximants 

In the previous section we calculated e expansions for the operator dimensions of the cubic 
Sp{N) theory in d = 6 — e. We may use the following Fade approximant to estimate the 


®This critical exponent is related to the derivative of the beta function in the 0{N) theory with quartic 
interaction; see 24 for more details on the comparison to the large N results. 
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Figure 3: Different Fade approximations of for the N = 2 model. 


behavior of operator dimensions as we continue d: 

Aq + Aid + A2d^ + ... + Amd^ 


PadeTT^^TT, [d] 


1 + Bid + ^2^2 + ... + Bnd"^ 


(2.7) 


For A = 2, by demanding that the e expansion have the same behavior as (2.4), we can £x 


four coefficients in the Fade approximant. Here we use Padei^ 2 , and Pade 2 ,i, to obtain the 
following estimate for the operator dimension in d = 5: 


A 


N=2 


= A 


N=2 


1.467 with Padei ^2 

1.459 with Fades,! • 


( 2 . 8 ) 


As expected, this is below the unitarity bound in d = 5, which is A = |. The plots of 
different Fade approximants are shown in Fig. For the next primary operator, whose 
scaling dimension is given in (3.11), using the Fades,! approximant we estimate A_|_ ^ 3.35 
in d = 5. In four dimensions we estimate A_|_ 2.7; this suggests that the N = 2 model 

does not have a free field description near d = 4. 

In the S'p(4) case, we may assume for 4 < d < 6 that the IR fixed point of the cubic 
theory (1.3) is equivalent to the UV fixed point of the quartic theory (1.1 )|^ If so, then not 
only do we know the 6 — e expansion of A^ to 0{e^), we can also use the quartic theory 


®This equivalence holds for Sp{N) models with sufficiently large N, but it is not completely clear if it 
applies for A = 4. 










result in (i = 4 + e, which is known to (9(e^): 

in (i = 6 — e 
in h = 4 + e 

(2.9) 

Together, these expansions allow us to determine ten coefficients of the Fade approximant. 
Disregarding the approximants that have a pole, we obtain the estimate 1.478 in 

d = 5. 


2 - 0.529827e - 0.0126197e2 + 0.0157244e3 

^x= { 

1 + 0.5e - 0.03125e2 + 0.015625e3 - 0.0730951e^ + 0.0195503e® 


2.2 Sphere free energies 


It is also interesting to compute the sphere free energy at the IR hxed point in d = 6 — e. 
The leading order term in the e expansion for the corresponding 0{N) models was computed 
. Sending N —>■ —N, we hnd the following result for F = sin(^) logZgd in the cubic 


m 


40 


Sp{N) models 


Fir = Fj 


uv 


TT {g*,Y-3N{glY 


-e + 0(e^), 


( 2 . 10 ) 


17280 (47r)3 

where -^uv = (1 — N)Fs, and Fg the value corresponding to a free conformal scalar. Plugging 
in the explicit solutions for the hxed point couplings gl, g^, it is straightforward to verify 
that, for all > 2, we have 

Fuv>Tir. (2.11) 


For instance, for N = 2 we hnd 


Fm — —Fg 


71 


43200 


+ 0{e^). 


( 2 . 12 ) 


Note that for the Fisher model (3^ of the Lee-Yang edge singularity, corresponding to Y = 0 


and imaginary g^, the inequality (2.11) does not hold. For general non-unitary theories the 
inequality does not have to hold; remarkably, it does hold for the Sp{N) models with positive 

N. 


Similarly, using the results of 40 for the quartic 0{N) theory in d = 4 — e, we can 


compute the sphere free energy at the IR hxed point of the model (1.1) in d = 4 — e. We 
hnd 

^ A'(lV- 2 ) , , „,4, 


Tquartic Tfree 


576 (Y-8)2 


+ Oie^) , 


(2.13) 


9 









where Ffree = —NFg. We observe that Ffree > -^quartic is satisfied for all iV > 2. 

The case = 8 is special and needs to be treated separately. Here the one-loop term in 
the beta fnnction vanishes, and we have 

A =-eA + + ©(A-). (2.14) 

Therefore, at the IR fixed point and we get 

-^quartic = -^free ~ . (2-15) 


3 Symmetry enhancement for N = 2 

Let ns write the cnbic Sp{2) model in terms of a real scalar a and a single complex anti- 
commnting fermion 6: 


S = j if X (d^Od^O + t + rjioOO + ^S2fj ■ 


(3.1) 


For g 2 = 2gi, i.e. the fixed point relation (2.3), this action possesses a fermionic symmetry 


with a complex anti-commnting scalar parameter a 


69 = aa , 69 = aa , 6a = —a9 + a9 


(3.2) 


As a conseqnence of this symmetry, the scaling dimensions of a and 9 are eqnal, as seen 


explicitly in eq. (2.4). This complex fermionic symmetry enhances the Sp{2) to OSp{l\2), 
which is the smallest snpergronp. The fnll set of snpergronp generators can be given in the 
form 


o"-5(4-4i 

■'•=•4 ■'-•'I- '■4(4-4i 


(3.3) 


lOpor e < 0 the interacting theory has a UV fixed point. So, in this case F^ee < Rjuartic, again in agreement 
with the conjectured F theorem. 
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and it is not hard to check that they satisfy the algebra of OSp{l\2): 

J^] = ±J^, [J+, J-] = 2J3 

[J“,Q±] =±iQ±, [J±,QT]=-Q± 


(3.4) 

(3.5) 

(3.6) 


We expect the operators of the theory to form representations of 05^(112). For example, 
let ns stndy the mixing of the Sp{2) singlet operators and 99. Setting = 2 in (2.5), we 
hnd the scaling dimensions of the two eigenstates: 

The hrst of these dimensions corresponds to the conformal primary operator 


0+ = a^ + 299, 


(3.8) 


which is invariant nnder all the 05p(l|2) generators. The second corresponds to 0_ = 
cT^ + 99, which is a conformal descendant becanse at the fixed point it is proportional to 
dfj,d^a by eqnations of motion. Indeed, we hnd A_ = + 2. 

Continnation of the resnlts for the cubic model with 0Sp{l\2) global symmetry to hnite 
e points to the existence of such interacting CFTs in integer dimensions below 6. In it 
was argued that the g —)■ 0 limit of the g-state Potts model is described by the 05^(112) 


sigma model 46,47. There is evidence that the upper critical dimension of the spanning 


forest model is 6 39 , and the 6 — e expansions of the critical indices are 39,48 


2 269 - 702C(3) 3 , 4, 

V = -e-e-e + 0{e ), 

' 15 225 16875 ^ ^ 


1 _ 7 
1 1 


3 30 27000 ^ ^ 


(3.9) 

(3.10) 


Remarkably, the operator dimensions we have calculated (2.4) and (3.7) agree with these 
expansions upon the standard identihcations 


d 


V 


o 1 + o ) 


A_i_ = d — u 


-1 


(3.11) 


Similarly, we can match the 6 — e expansions of the sphere free energies. For the g-state 
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Potts model it is not hard to show that 


Fq — {q~ ^)Fs + 


7r(g-l)(g-2) 2 


8640(3g - 10) 


6 ^ + 0 ( 6 ^), 


(3.12) 


and for q = 0 this matches the F of our 05^(112) model, (2.12). These results provide strong 


evidence that the 0Sp{l\2) symmetric IR fixed point of the cubic theory (3.1) describes the 
g —)■ 0 limit of the g-state Potts model. 


Numerical simulations of the spanning-forest model 39 , which is equivalent to the g —)■ 0 


limit of the ferromagnetic g-state Potts model, indicate the existence of second-order phase 


transitions in dimensions 3, 4 and 5. The critical exponents found in 39 are in good 


agreement with the Pade extrapolations of 6 — e expansions exhibited in section 2.1. This 
provides additional evidence for the existence of the critical theories with OSp{l\2) symmetry. 
It would be of further interest to find the critical statistical models that are described by the 
Sp{N) invariant theories with N > 2. 
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